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ABSTRACT 

Context. The Galactic Center offers unique opportunities to study stellar and bow-shock polarization effects in a dusty environment. 
Aims. The goals of this work are to provide near-infrared (NIR) polarimetry of the stellar sources in the central parsec at the resolution 
of an 8m telescope for the first time, along with new insights into the nature of the known bright bow-shock sources. 
Methods. We use adaptive-optics assisted observations obtained at the ESO VLT in the H- and Ks-band, applying both high-precision 
photometric methods specifically developed for crowded fields and a newly established polarimetric calibration for NACO to produce 
polarization maps of the central 3"xl9", in addition to spatially resolved polarimetry and a flux variability analysis on the extended 
sources in this region. 

Results. We find foreground polarization mainly parallel to the Galactic plane, with average values of (4.6 ± 0.6)% at 26° ± 6° (Ks- 
band) and (9.3 ± 1.3)% at 20° ± 6° (H-band) in the center of the field-of-view (FOV). Further away from the center, we find higher 
polarization degrees and steeper polarization angles: (7.5 ± 1.0)% at 11° ± 6° (Ks-band) and (12.1 ± 2.1)% at 13° ± 6° (H-band). 
Ph/Pks peaks at 1.9 ±0.4, corresponding to a power law index for the wavelength dependency of o- = 2.4 ±0.7. These values also vary 
over the FOV, with higher values in the center This is indicative of the influence of local effects on the total polarization, possibly 
dichroic extinction by Northern Arm dust. The two extended sources IRS 21 and IW show similar intrinsic polarization degrees of 
6.1 resp. 7.8% (Ks) and 6.9 (H, only IW) at polarization angles coincident with previous NIR and mid-infrared (MIR) findings, both 
in total and spatially resolved. The spatial polarization pattern of both sources points to scattering on aligned elongated dust grains as 
the major source of intrinsic polarization, and matches the known orientation of the magnetic field. Our data also allow us to separate 
the bow shock of IRS 21 from the central source for the first time in the Ks-band, finding the apex north of the central source and 
determining a standoff distance of ~400 AU, which matches previous estimates. This source also shows a ~ 50% increase in flux in 
the NIR over several years. 
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(^ 1. Introduction Milky Way spiral arms along the line-of-sight (LOS). The 11.5 

T— I fim polarization, however, was found to be almost perpendicu- 

'— I ' The center o f the Milk y Way is located at a distance of ~ 8.0 kpc lar to the galactic plane and was therefore classified as intrinsic. 

^ ; dGhez et all [2008; Gill essen et"dl l2009l) . It contains a nuclear The values found for the L-band showed intermediate values, 

stellar cluster (NSC) with a ~ 4.0x 10^ Mp super -massive black which was attributed to a superposition of both effects. The GC 
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^ hole at its dynamical center jlEckart et all 120 021: [Schodel et al.. therefore offers the possibiUty of studying both interstellar po- 



H i I2OO2I l2003t iGhez et all [20031 120081: iGillessen et al.Ll2009l) ."fhrs larization and th e prope rties of intrinsically polarized sources. 

5t , cluster shows similar properties as the NSCs found at the dy- iKobavashi et alj (Il980l) conducted a survey of K-band polar- 

namical and pho tometric centers of other galaxies (Boker, 20l2 ization in a much wider field of view (7'x7'), finding largely 

ISch6delLl2010cl) . uniform polarization along the galactic plane. iLebofskv et al.l 

Less than a decade after the first near-inf rared (NIR) imaging ob- (1982) confirmed these findings in the H- and K-band for 17 

servations of the Galactic center (GC) bv lBeckUn & NeugebaueJ sources in the central c l uster. The latest H-band survey was con- 

(I1968'). the polarization of 10 sources within <2 pc of Sagittarius ducted bv iBailev et~an (Il984j) . who examined -10 sources with 

A* wa s measured by Capps & Knacke ( 1976); Knacke & Capps a 4.5" resolution, finding similar results. While these surveys 

(Il977h . observing in the K-, L- and 1 1.5 yum-band. These obser- could only resolve a small number of sources in th e central re- 

vations revealed similar polarization degrees and angles for the gion, h igher resolution observations (0.25") enabled .Eckart et al.1 

four sources observed in the K-band, with the polarization angles ( Il995h to meas ure the polarizati on of 160 sources in the central 

roughly parallel to the galactic plane (-4% at 15° East-of-North, 13"xl3", while lOttet all d 19991) examined ~ 40 bright sources 

while the projection of the galactic p lane is at a positio n angle of in the centtal 20"x20" at 0.5" resolution. These two studies 

~3 1 .4° at the location of the GC, see iReid et aUl2004 . This was confirmed the known largely uniform foreground polarization, 

interpreted as polarization induced by aligned dust grains in the but already revealed a more complex picture: individual sources 
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Fig. 1. Ks-band field-of-view, one channel of the Wollaston prism. 



showed different polarization parameters, such as a significantly 
higher polarization degree (IRS 21). 

These results akeady illustrated that in addition to the fore- 
ground polarization, intrinsic polarization plays a role in the Ks- 
band. Especially sources embedded in the Northern Arm and 
other bow-shock sources (see Fig|2]for an overview of the bright 
bow-shock sources in the central 20") show signs of intrinsic 
polarization. Among these objects, IRS 21 shows the strongest 
total Ks-ban d polarization of a bright GC source detected to 
date (Eckart et al.. 119951: lOtt et al.L I1999l ~ 10-16% at 16°). 
iTanner et al.. (2002) described IRS 2 1 as a bow-shock most likely 
created by a mass-losing Wolf-Rayet star. The observed polar- 
ization is a superposition of foreground polarization and source 
intrinsic polarization. It is still unclear what process(es) cause 
the latter component: (Mie)-scattering in the dusty environment 
of the Northern Arm and/or emission from magnet ically aligned 
dust. The latter is known to occur at 12.5 /im dAitken et al.L 
Il998h . but should be negligible at shorter wavelengths. 
In general, three effects can produce polarized NIR radiation 
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Fig. 2. Lp-band image of the innermost 20" of the GC (ESQ VLT 
NACO image, program 179.B-0261(A)). Prototypical (bow-shock) 
sources and location of the Northern Arm are indicated. Stellar sources 
are marked with stars, bow-shocks are indicated by triangles. 



in GC sources: (re)emission by heated, non-spherical grains, 
scattering (on spherical and/or aligned non-spherical grains) and 
dichroic extinction by aligned dust grains. The first two cases can 
be regarded as intrinsic to the source for our purposes, thereby 
allowing conclusions about the source itself and its immediate 
environment, while the third effect is the result of grain align- 
ment averaged along the LOS. For sources enclosed in an op- 
tically thick dust shell, dichroic extinction can also contribute 
significantly as a local effect (see e.g. Whitney & Wolff , 2002). 
As the basic mechanis m that could cause t h e obse rved large- 
scale grain alignment, iDavis & GreensteinI ( Il951h suggested 
paramagnetic dissipation, which basically aligns the angular mo- 
mentum of spinning grains with the magnetic field. But even 
almost 60 years later, the problem of grain alignment is by no 
means completely solved, and it remains difficult to reach exact 
conclusions for dust parameters and magnetic field strength, but 
at least determining the magnetic field orientation is possible. If 
the parameters change along the LOS, this further complicates 
the issue. See e.g. Purcell et al. (1971), Lazarian (2003), cite- 
lazarian2007 and references therein for a review of the different 
possible causes of grain alignment expected to be relevant in dif- 
ferent environments. 

This makes it possible to use polarimetric measurements to 
map at least the direction of the magnetic fields responsible 
for dust alignment through the Davis-Greenstein effect, as e.g. 
iNishivama et al.l (l2009l 1201 Ol) showed for the innermost 20' 
resp. 2° of the GC, but these studies did not cover the central 
parsec owing to insufficient resolution. 

Observations of the Galactic Center suffer from strong extinc- 
tion caused by dust grains on the LOS, with values of up to 
Ay - 40 mag at optical wavelengths (or even higher values 
of up to 50 mag if a steeper extinction law is assumed) and 
still around 3 mag in the Ks-band (e.g. Scoville et al.L 120031 : 
ISchodel et al.L l2010bh . The aligned interstellar dust grains re- 
sponsible for the polarization cause extinction as well, but non- 
aligned grains can also contribute. Therefore, the same particles 
are no t necessarily responsible for both effects (e.g. Martin et al], 
Il990h . Universal po wer laws have bee n claimed for NIR e xtinc- 
tion bv lDraind(IT989.) and polaiization (iMartin et al.LI 19901) . who 
also showed that the law applicable to polarization in the opti- 
cal domain dSerkowski et all Il975h is a poor approximation in 
the NIR. The power law indices presented in these works for the 
extinction and the polarization power law are almost the same 
(1.5-2.0). It also appears that there is a correlation between the 
measured extinction of a n intrinsically unpolar ized source and 
foreground polarization dSerkowski et all 11975 ). but this rela- 
tion is quite complex. In the Ught of new results for the extinc- 
tion law, which seem to deviate consistently from the Draine law 
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Fig. 4. H- (green) and Ks-band (red) lightcurve of the known non- 
variable source IRS 16C. 
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Fig. 3. Total relative flux errors of the sources detected in the Ks- (up- 
per frame) resp. H-band (lower frame). Only the 0° channel is shown 
here, but the error distribution is very similar in the other three chan- 
nels. The values are based on variations between the fluxes measured 
for the same source at different dither positions. The red lines denote 
the 16/18 mag resp. 3/6% upper limits (see Appendix IaTi. 



(e.g. lGosling et all 120091: FSchodel et a l.', '2010b'), new polarimet- 
ric measurements may be useful to further clarify the relation 
between extinction and polarization. The central parsec of the 
GC is a well-suited but challenging environment to study this re- 
lation, because it contains a huge number of sources that exhibit 
large variances in exti nction (1-2 mag, see lBuchholz et al.Ll2009t 
ISchodel et all 1201 0b'), which is produced along a long line-of- 
sight by a great number of dust clouds with possibly different 
grain alignment and composition. We are using a new Ks-band 
extinction map of the central parsec recently been presented by 
ISchodel et all (I2010bh . 

The aims of this work are to present the first map of stellar H- 
and Ks-band polarization of the central few arcseconds, a highly 
crowded region with more than 10 sources per arcsec^. In order 
to avoid crowding and to resolve the sources individually, the 
resolution of an 8m class telescope is required. Especially in the 
H-band, no polarization measurements exist so far for individ- 
ual resolved sources in the central parsec. In addition, we will 
present the first spatially resolved polarimetric measurements on 
IRS IW and IRS 21, also resolving the bow-shock structure of 



IRS 21 for the first time in the Ks-band. Furthermore, we will 
conduct a variability analysis on the extended sources IRS IW, 
5, 10 and 21 in the H-, Ks- and L-band. 

In ^ we discuss the photometric methods used and the cali- 
bration applied to the data. Our results will be presented in ^ 
followed by a summary and discussion of their implications in 



2. Observation and data reduction 

2.1. Observation 

The polarimetric datasets used here were obtained using the 
NAOS-CONICA (NACO) instrument at the ESQ VLT unit tele- 
scope 4 on Paranal in June 2004 (H-band broadband filter, pro- 
gram 073.B-0084A, dataset 1, see TabO]) and May 2009 (Ks- 
band broadband filter, program 083.B-0031A, dataset 2). We 
also used several additional polarimetric Ks-band datasets con- 
tained in the ESO archive (datasets 3-15l3that covered IRS 21 as 
well through a rotated field-of-view (FOV). Please see Tab[3]for 
details on the individual observations. For our variability study 
of the bow-shock sources we used NACO imaging data taken 
between June 2002 and May 2008 (H-, Ks- and L'-band data 
contained in the ESO archive). 

The seeing during the polarimetric Ks-band observations was 
excellent with a value of ~0.5", while conditions during the H- 
band observations were less optimal with a seeing of ~0.8". We 
used the bright super-giant IRS 7 located about 6" north of Sgr 
A* to close the feedback loop of the adaptive optics (AO) sys- 
tem, thus making use of the infrared wavefront sensor installed 
with NAOS. The sky background was determined by taking sev- 
eral dithered exposures of a region largely devoid of stars, a 
dark cloud 713" west and 400" north of Sgr A*. The Wollaston 
prism available with NACO in combination with a rotatable half- 
wave plate was used for the polarization measurements. The 
two channels produced by the Wollaston prism (0° and 90°), 
combined with two orientations of the half-wave plate (0° and 
22.5°), yielded four sub-images for each of several dither posi- 
tions along the east-west axis. In total, we were able to cover 
a field-of-view of 3"xl6.5" (Ks-band) respectively 3"xl9" (H- 
band), corresponding to 0.12 pc x 0.66 pc and 0.12 pc x 0.72 
pc, respectively (see FiglTJ- 



' Based on observations collected at the European Organization for 
Astronomical Research in the Southern Hemisphere, Chile 
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Fig. 5. Ks-band polarization map of stars in the Galactic center. Only reliably measured polarization values are shown here. Blue circles and green 
lines: dataset 2. Black circles and red lines: dataset 4. The diameter of the circles corresponds to the brightness of the source. For comparison, we 
also show preliminary results of new observations taken in March 201 1 (denoted by red circles and black lines) . Typical errors are indicated by 
error cone (lower center). Thin black lines in the background denote magnetic fields determined from MIR data dAitken et allll99§ . , based on a 
1.5" beam). The brightest sources are also indicated. 



All images were corrected for dead/hot pixels, sky-subtracted 
and flat-fielded. It is essential for a correct calibration that the 
flat-field observations are taken through the Wollaston, because 
the flat-field shows variations caused by different transmissivity 
of flie channel s as well as efifects o f the inclined mirrors behind 
the prism (see IWitzel et al.L 12011 . for details on these effects). 
Using a flat-field taken without the Wollaston prism can produce 
offsets in the order of 2% and 10° in the measured polarization 
parameters. FiglT] shows the region covered by the final mosaic 



that was used for the polarimetry. 

The data used here were originall y take r i to examine flares 
of Sgr A* (see eg. lEckart et'aLl l2006t iMever et al.L 120061: 
IZamaninasab et al1.l2010l) . and this purpose requires the highest 
possible time resolution. This is the reason why differential po- 
larimetry could not be applied. This technique can eliminate or 
reduce many instrumental effects and thus increase precision by 
rotating the imager by 90° between exposures and thus canceling 
these effects out (effectively switching the and 90° channels). 
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Fig. 6. Ks-band polarization degrees (plotted on logarithmic scale, upper left) and angles (lower left) of stars in the Galactic Center (2009-05-18, 
dataset 2). The red line denotes the fit with one Gaussian distribution, while the green resp. blue lines denote the fit with a double Gaussian (green: 
sum, blue: individual Gaussians). Upper right: relative errors of the polarization degrees. Lower right: absolute errors of the polarization angles. 



But this also means more time is needed for each exposure, be- 
cause two images have to be taken for one data-point and the 
actual rotating takes time as well. 



2.2. Photometry 

2.2.1. Deconvolution-assisted large-scale photometry 

For all datasets, the individual exposures were combined to a 
mosaic. All photometry was conducted on these mosaics. 
Accurate photometry is crucial for polarimetry, especially when 
the polarization of the targets is of just a few percent. This is the 
case for the sources in the central parsec in the H- and Ks-band, 
and the effects of crowding and variations of the point-spread- 
function (PSF) over the FOV complicate photometry even more. 
For the very bright sources in the FOV, such as the IRS- 16 and 
IRS-1 sources and the extended sources IRS-21 and IRS-IW, 
crowding is not a problem, but saturation can lead to additional 
complications. These effects have to be countered effectively in 
order to achieve low photometric eiTors. 

We therefore adopted a photometric method recently presented 
by ISchodell (201()ah : first, we used the StarFinder IDL code 



dPiolaiti et al. 



20001) to repair the cores of saturated sources 



(only necessary for some very bright sources in the H-band im- 
age) and extract a PSF on the full image from sufficiently bright 
and isolated sources. For this first step, the most suitable source 
would be the guide star IRS 7 itself, since it is several magni- 
tudes brighter than any source within several arcseconds, but 



this source was not covered by the FOV of our dataset. We used 
a PSF determined from several IRS 16 and IRS 1 sources in- 
stead. Unfortunately, the most suitable of these sources were 
only contained in the FOV of the Ks-band data, but not cov- 
ered in the H-band. Since this process is d esigned to determine 
the faint wings of the PSF accurately (see Buchholz et al.L 120091: 
ISchodel , 2010a), it works worse the fainter the initially used 
PSF sources are. 

We used this PSF for a linear deconvolution, i.e. a division in 
Fourier space, followed by the application of a Wiener filter We 
then applied local PSF fitting photometry by using StarFinder on 
overlapping sub-frames of th e deconvo l ved im age (~3"xl0.6", 
using wider sub-fields than ISchodel I (1201 Oah . this generates 
large overlapping regions and ensures that enough bright sources 
are contained in each subfield for the PSF estimation). This 
method significantly reduces source confusion and also mini- 
mizes systematic errors introduced by variations of the shape of 
the PSF over the FOV due to anisoplanasy. Despite the relatively 
narrow FOV, this effect still occurs because the distance and the 
angle to the guide star change considerably over the field. While 
using naiTower sub-fields might counteract this, the availability 
of sufficiently bright PSF stars takes precedence, since this fac- 
tor is the primary limit for the quality of the photometry. 
The resulting fluxes in all four channels were normalized to an 
average of one over all channels for each source and then merged 
to a common list of sources detected in all channels of each sub- 
frame. The sub-frame lists were then merged to a common list 
of all detected sources. 
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Fig. 7. Ks-band polarization degrees (plotted on logarithmic scale, upper left) and angles (lower left) of stars in the Galactic Center (2007-04-03, 
dataset 4). The red lines denotes fits with a Gaussian distribution. Upper right: relative errors of the polarization degrees. Lower right: absolute 
errors of the polarization angles. 



Fig[3]show the photometric uncertainties in the H- and Ks-band. 
We only used sources brighter than 16 mag (Ks-band) respec- 
tively 18 mag (H-band) and with relative photometric errors of 
less than 3% resp. 6% in the subsequent analysis. The lower 
brightness limit was chosen to avoid problems with insufficient 
completeness and unreliable photometry, because the errors in- 
crease drastically for fainter sources. Please see Appendix |A] 
for details on our error estimation and the reasons for the error 
threshold. 



2.2.2. Photometry on extended sources 

Applying the PSF fitting algorithm to extended sources leaves 
large residua, and a simple core-subtraction with a stellar PSF 
does not counter any distortions of the extended component pro- 
duced by the atmosphere and the telescope. But high precision is 
necessary here as well. We therefore used a different approach: 
we deconvolved a small region of the mosaic images with the 
Lucy-Richardson algorithm (the "ringing" produced by the lin- 
ear deconvolution complicates the photometry of extended fea- 
tures), using the PSF estimated from the complete mosaic. This 
resulted in a much clearer view of the extended features (see 
FigfTSll. while other sources in the vicinity appear point-like. 
We shifted the resulting images to a common reference frame. 
We then covered the extended features with overlapping aper- 
tures (~ 27 mas radius), measuring the flux in each aperture. 
From these fluxes we determined the polarization at the position 
of the apertures. For comparison, the same method was applied 



to the PSF used for the deconvolution (see Appendix iBli. In ad- 
dition, we determined the total polarization of the two known 
extended sources in our FOV, IRS IW and IRS 21, by covering 
them with apertures of 0.25" radius. This allowed a comparison 
to previous observations. 

In order to study a possible flux variability of the extended 
sources (see ^3.6l l, we used 45 Ks-band and 38 H-band NACO 
datasets contained in the ESO archive (taken between 2003 and 
2008). Where the bright point-sources were saturated, we used 
StarFinder to repair their cores. We then applied aperture pho- 
tometry to the exten ded sourc es and sever al known bright non- 
variable sources (see lOtt et an.[T999.;Rafelski et aI.Li2006.) . with 
additional apertures placed on regions devoid of stars close to the 
targets to obtain a background estimate. This background was 
subtracted from the recovered fluxes and the result normalized 
to the total flux of the chosen non-variable calibration sources 
for each dataset (the IRS 16 sources except 16SW and 16NE 
and the IRS 33 sources; IRS 16NE is strongly saturated, while 
IRS 16SW has been described as an eclipsing binary and there- 
fore shows strong fl ux variabilities, making it unsuitable as a 
reference source, see lOtt et aI.llI999HRafeIski et aI.Ll2006h . We 
estimated the uncertainty of each flux measurement by repeating 
the aperture photometry on the individual non-mosaiced images 
of each dataset (which can be assumed to be independent mea- 
surements), and adopted the standard deviation of the recovered 
fluxes as the total flux error. Especially for bright and isolated 
sources, this yields very small statistical errors on the order of 
less than 0.5%. We still find systematic variations of the mea- 
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Fig. 8. Comparison between polarization degrees (left frame) and angles (right frame) measured on sources common to dataset 2 (2009) and 
dataset 4 (2007). Plotted are histograms of — resp. O2 - 64, with Gaussians fitted to both values (red lines). 



sured fluxes on the order of about 5% between the epochs even 
for sources known to be non-variable on the timescale used here 
(see FigHJi. The latter value provides a more realistic estimate of 
the photometric accuracy, so we introduced an additional flux er- 
ror of 5% for all data-points to be able to separate real variability 
from noise. 



2.3. Polarimetry 

We determined the polarization degree and angle of each source 
by converting the measured normalized fluxes into normalized 
Stokes parameters: 



/ = 1 

g _ /o - /90 
/o + /90 

jj _ /45 - /l35 
/45 + /l35 

y = 0. 



(1) 

(2) 



Because NACO is not equipped with a j plate, it was not pos- 
sible t o measure circular polarization. However, iBailev et al.l 
(1198 4*) showed that the circular polarization of sources in the 
GC is at best very small, so we assume here that it can be ne- 
glected and set to at our level of accuracy. Polarization degree 
and angle can then be determined as 



P = Vfi^ + u^ 
6* = 0.5 X atari — 



(3) 
(4) 



Errors for I, Q, U and subsequently p and 9 were determined 
from the flux errors. In order to determine whether or not the 
polarization of a source was determined reliably, we calculated 
the normalized fluxes that would be expected for the determined 
values of p and 6. The difference between these values and the 
measured fluxes was then compared to the photometric errors 
of each data-point. The source was only classified as reliable if 
the root-mean-square of the deviations did not exceed the root- 
mean-square of the relative photometric errors. The subsequent 
analysis is only based on these higher quality sources. 



2.4. Calibration of tlie measured polarization 

A first comparison of our measurements to known values 
(Knacke & CaDDS.'1977t [Bailev et allfTosilEckart et"ailll995l: 
Ott et al., 1999) revealed significant offsets, especially in the 
polarization angles. Instead of the expected orientation along the 
galactic plane (oriented ~3 1 .4° East-of-North, in the following, 
positive angles should be read as East-of-North, negative as 
West-of-North), we found orientations of the polarization vec- 
tors of about -5°. The problem that surfaces here is that NACO 
was not specifically built for polarimetry, so instrumental effects 
like this can be expected and have to be countered by a special 
calibration. Also, the observation technique was developed for 
the study of short-term variabilities (flares) of Sgr A* and not 
for high-precision polarime try on stellar sources. 
To reduce this problem, IWitzel et alJ (1201 ih developed an 
analytical model of the behavior of polarized light within 
NACO. It consists of Miiller matrices to be applied to the 
measured Stokes vector of each source, which then yields the 
actual Stokes parameters of the source. In general, any optical 
effect such as reflection, transmission, polarization etc. can 
be described by a Miiller matrix. A combination of effects as 
it occurs here is then represented by a multiplication of the 
individual Miiller matrices. If the necessary material constants 
and the construction of the instrument are known, the resulting 
matrix can be us ed to significantly r educe systematic offsets 
and uncertainties. IWitzel et alJ (1201 ih show that by using this 
method, the systematic uncertainties of polarization degrees and 
angles can be reduced to ~1% and ~5°. Applied to our data, 
this limits the final accuracy by the photometric uncertainties 
instead of by instrumental effects. For a detailed d escription of 
the mo del and the Miiller matrices themselves, see IWitzel et al.l 
(1201 Ih . 

Utilizing this new calibration model means that an actual direct 
calibration can be ac hieved for the first t ime at this resolution. 
Previous studies like lEckart et alJ (Il995h and Ott et al. (199i) 
had to adopt a ca l ibratio n based on reference values taken from 
lKnacke&Cappsl(ll977h . 



2.5. Correcting for foreground polarization 

It can be assumed that the total effect of the foreground polar- 
ization can be treated as a simple linear polarizer with a certain 
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Fig. 9. H-band polarization map of stars in the Galactic Center (dataset 1). Only reliably measured polarization values are shown here. Typical 
errors plotted as error cone (upper right), adopted from Ks-band data. 



orientation 6fg and efficiency p/g. This can be described by a 
Miiller matrix: 



* obs 
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with Sobs as the observed total Stokes vector and 5,,,, as the 
Stokes vector of the intrinsic polarization. Munip) is the Miiller 
matrix describing a linear polarizer, producing a maximum of 
polarization along the North-South-axis: 



Mlinip) = 



. This matrix has to be rotated to the appropriate angle by multi- 
plying it with M,-oriO'), a standard 4x4 rotation matrix: 
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3. Results and discussion 

3.1. Ks-band polarization 

2009 data (dataset 2) We were able to measure reliable polar- 
ization parameters for 194 sources brighter than 16 mag in our 
main 2009 Ks-band dataset. For fainter sources, the photometric 
uncertainty becomes too large to determine the polarization reli- 
ably, and this limit also helps to avoid problems with insufficient 
completeness and source crowding, which could generate a bias 
in the averaged values. The subsequent analysis and the compar- 
ison to the H-band data is only based on this dataset. 
The polarization angles in the central arcseconds mostly follow 
the orienta tion of the galact ic plane within the uncertainty lim- 
its (31.4° (lReidetal.Ll2004l) . while we find angles of -25-30° 
for our sources). Toward the eastern edge of the FOV, we find 
slightly steeper angles (~5-15°, see Fig|5]l. A few sources west 
of Sgr A* also show similar steep angles, but there are too few 
reliable sources there to allow any conclusions. Here and in the 
following, we use the term "steeper angle" to denote angles with 
an absolute value closer to 0°, so that an angle of 10° would be 
steeper than one of 30°. The polarization angle can vary between 
-90° and 90°, in the way that 91° (East-of-North) corresponds to 
-89°. 

The distribution of the polarization angles (see Fig|6] lower left 
frame) can be fitted with a single Gaussian, peaked at 20° with 
a FWHM of 30°. Using the FWHM as a measure for the un- 



certainty (with (T = 



FWHM 



), this yields = 20° + 13° 



. 0' - 90° + 0fg has to be used in the rotation matrix, because we 
define the polarization angle as the angle where we measure 
the flux maximum, while the angle of reference for the Miiller 
matrix describing the linear polarizer is the angle where the max- 
imum in absorption occurs. 

For each source to which the depolarization matrix was applied, 
we used the average of the polarization parameters of the sur- 
rounding point sources as an estimate for 0fg and pfg. 
The resulting matrix can then be inverted and multiplied with the 
calibrated observed Stokes vector of a source to remove the fore- 
ground polarization and leave only the intrinsic polarization. We 
applied this method to the total polarization and the polarization 
maps of the extended sources (IRS 21 and IRS IW), to isolate 
their intrinsic polarization pattern. The relevance of the results 
depends heavily on the accuracy of the foreground polarization 
estimate. 



2 V2x/n(2) ^ 

Using a fitting function with two Gaussian peaks yields a sig- 
nificantly lower ;i'^ (by a factor of 3). The two peaks are fitted at 
10° ±4° (FWHM of 10°) respectively 6I2 = 26°+8° (FWHM 



I 
of 20°). Considering the uncertainties of the polarization angles 
shown in the lower right frame of Fig|6] which are on the order 
of up to 15°, it can be questioned if these two peaks are indeed a 
real feature, with the distance between the peaks on the order of 
these typical eiTors. 

The polarization degree also appears to vary over the field, with 
values of 4-5% in the central region and 8-10% toward the east- 
ern edge (and for some western sources, but with the same caveat 
as for the polarization angle). Especially sources in the area 
around the IRS 1 sources show these higher polarization de- 
grees. We fitted the logarithms of the polarization degrees with 
a Gaussian (peaked at (5.1 + 1.7)%, FWHM of 4.0%, see Fig|6l 
upper left frame), and like the single Gaussian fit to the polar- 
ization angles, the fit was quite poor Repeating the fit with a 
double Gaussian yielded two peaks at (4.6 + 0.8)% respectively 
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Fig. 10. H-band polarization degrees (plotted on logarithmic scale, upper left) and angles (lower left) of stars in the Galactic Center (dataset 1). 
The red line denotes the fit with one Gaussian distribution, while green resp. blue lines denote the fit with a double Gaussian (green: sum, blue: 
individual Gaussians). Upper right: relative errors of the polarization degrees. Lower right: absolute errors of the polarization angles. 



(7.7 + 1.2)%, with FWHMs of 1.8% resp. 2.8% and a signifi- 
cantly better ;t'^ (by a factor of 8). The relative uncertainties of 
the polarization degree are on the order of up to 30% (see Fig|6] 
upper right frame), and this limits the confidence in the two fit- 
ted peaks. 

Comparing the two fitted Gaussian distributions for both param- 
eters, we find a similar number of stars contained in the 10° and 
the 7.7% peak (-25-30%), resp. the 26° and 4.6% peak (-70- 
75%). This confirms the general trends found in Fig |5] and indi- 
cates that the fitted peaks indeed correspond to a real feature. 
We plotted the Ks-band polarization degrees versus the polariza- 
tion angles in FigHT] left frame, and there appears to be a trend 
that higher polarization degrees coincide with steeper polariza- 
tion angles, despite the large errors. 

2007 data, rotated FOV (dataset 4) We measured reliable 
polarization parameters for 186 sources brighter than 16 mag. 
This dataset has a lower Strehl ratio (on average) than the 
2009 data (22% compared to 27%). We used this dataset only 
for comparison with the main Ks-band data (dataset 2, with 
the main purpose to determine the presence of unaccounted 
instrumental effects over the FOV), since no H-band data with 
this FOV is available. 

The trends we find here are similar to those in the main dataset: 
the polarization angles appear to be aligned with the Galactic 
plane (see Fig|5]l, but we do not observe the same shift in polar- 
ization angle toward the east of the FOV. We do, however, find 
an increase in polarization degree toward the southeast, similar 



to the increase found in the 2009 data toward the east. Both the 
distribution of the polarization degrees (logarithmic) and the 
polarization angles can be fitted with a single Gaussian with 
sufficient accuracy (see FigQ, with peaks fitted at (4.6 ± 2.1)% 
respectively 26° ± 8° (FWHM of 5% resp. 20°). The FWHM of 
the distribution of the polarization degrees is comparable to that 
found for the 2009 data, but a single Gaussian provides a much 
better fit here. The larger uncertainties would however lead to a 
bluiTing of the two Gaussians, if indeed two were present. 
The relative eiTors of the polarization degree (see Fig|6] upper 
right frame) mostly stay below 30%, with some outlier values 
of up to 40-50%. This exceeds the errors found for dataset 
2, but this can be expected because of the lower data quality. 
The errors found for the polarization angles are also larger on 
average than those measured for dataset 2. 
The sources in the overlapping area of dataset 2 and 4 for which 
reliable measurements could be obtained in both datasets show 
very similar polarization parameters (see Fig|8]l. We find that the 
polarization degrees and angles of 82% of the common sources 
agree within 1 sigma, viewing the parameters individually. 
Both parameters agree within 1 sigma for 69% of the common 
sources. For 93% of the common sources both parameters agree 
within 2 sigma, and 98% (all but one source) show an agreement 
within 3 sigma. 

We fitted the relation between the polarization degrees resp. 
the difference of the polarization angles of the sources in the 
common region of both datasets with a Gaussian (see Fig|8]l. 



We fitted the peaks at 
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1.2 ± 0.2 (FWHM of 0.5) and 
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Fig. 11. Polarization angle vs. polarization degree, plotted as logarithmic smoothed point density, with typical error indicated by red bars (upper 
right corners). Contours were plotted to guide the eye. Left frame: Ks-band, right frame: H-band 



02-04 ^ r ± 7° (FWHM of 17°). The uncertainties of these 
values provide an estimate for the general accuracy of the 
Ks-band measurements, so we plotted them as an error cone in 
FigS 

This reinforces the confidence in the measured values in 
both datasets, while it also provides a confidence limit. A 
3cr-deviation would correspond to a (relative) 60% offset in 
polarization degree and a difference of 21° in polarization angle. 
Unfortunately, the overlapping part of the FOVs is only about 
8 arcsec^, and no data with a different FOV exist for the 
regions that show an excess in polarization degree/angle in both 
datasets. 

FigH also shows the preliminary results of new polarization 
measurements (program 086.C-0049(A)) of a region north of 
the FOV of dataset 2. These measurements seem to link up well 
with the findings in dataset 2 and 4. A detailed analysis of these 
observations will be the subject of a future study. 



3.2. H-band polarization 

Reliable results could be obtained for 163 sources brighter than 
18 mag. The limit of 18 mag was chosen because this corre- 
sponds to the limit of 16 mag in the Ks-band, assuming a typi- 
cal H-Ks of ~2 (with H-Ks,„,„>,s,c ~0 and A^-Aks ~2, see e.g. 
ISchodel et all l2010bl) . The lower number of sources with reli- 
able polarization compared to the Ks-band can be attributed to 
the significantly lower Strehl ratio of the H-band data (0. 17 com- 
pared to 0.27 in the Ks-band) and the slightly different FOV, but 
the latter is a minor effect. We find polarization angles very sim- 
ilar to those in the Ks-band, but with a more uniform distribution 
over the FOV (see Fig|9]l- The distribution of the polarization an- 
gles can be fitted well with a single Gaussian, peaking at 20° ±8° 
(FWHM of 19°, see FigQUl lower left frame). Fitting this distri- 
bution with a double Gaussian produces a slightly better ;(f", but 
this can be expected for an increase in the number of fitting pa- 
rameters. A single Gaussian fits the distribution with sufficient 
accuracy, compared to the poor fit with a single Gaussian func- 
tion for the Ks-band polarization angles. We find typical errors 
of the polarization angle of up to 12° (see FiglTO] lower right 
frame). 

The polarization degree also appears to be quite uniform over 
the FOV, with typical values of 8-12%. Fitting the logarithms of 



the polarization degrees with a single Gaussian leads to a peak 
at (9.8 ± 0.7)% (FWHM of 1.7%), satisfyingly matching the 
data (see Fig[TOl upper left frame). We also fitted the data with 
two Gaussian peaks for comparison, but this only marginally im- 
proves the fit. We find relative uncertainties of the polarization 
degree on the order of up to 40% (see Fig[TOl upper right frame). 
Unfortunately, no high-resolution H-band measurements with a 
different FOV are available for a comparison. 
As for the Ks-band, we plotted the H-band polarization de- 
grees versus the polarization angles (see Fig[TT] right frame). No 
large-scale trend is visible, both parameters appear concentrated 
around -20° resp. -9%. 



3.3. Comparison to previous results 

We c ompared the polariz ation degrees and ang les of 30 sources 
from lEckartet all (11995 !) and 13 sources from lOtt et all (Il999h 
to the values we determined from our data. These two studies 
were calibrated based o n po larization parameters d etermined by 
iKnacke & CappsI (Il977h and lLebofskv et al.1 (Il982h . We find that 
87% of the Eckart sources agree with our values within 3 sigma 
in polarization degree (and 83% in polarization angle). For the 
Ott sources (with a slightly different FOV) we find that 77% 
agree with our own values within 3 sigma for both polarization 
degree and angle. Differences between these older studies and 
our own measurements can probably be attributed to the lower 
spatial resolution of the former. Both older studies show average 
polarization angles generally parallel to the Galactic plane, on 
average at 25° resp. 30°. We find similar results for our much 
larger number of sources. The average polarization degree is 
slightly higher, but this can be expected because the inclusion 
of IRS 7 with its polarization degree of only 3.6% in both older 
surveys considerably lowered the flux-weighted average that was 
calculated there. 

One other thing has to be considered here: we have for the first 
time applied an absolute polarimetric calibration to high an- 
gular resolution data. That we find such good agreement with 
data w here a relative calibration based on the .Knacke & CappsI 
(Il977h values was used increases the confidence in both this 
study and the results of the mentioned previous works. 
T he only direct possibility to compare our study with the results 
of iKnacke & CappsI (1 19771) exists in the vicinity of IRS 1 . There, 
we find a flux-weighted average polarization of 4.3 + 0.5% at 
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Fig. 12. Left frame: difference between H- and Ks-band polarization angle. Right frame: relation of H- to Ks-band polarization degree (logarithmic 
plot). The red line represents a Gaussian fitted to the histograms. 



14° ± 10° in a 3.5" aperture, which matches the value given in 
the older study (3.5 + 0.5% at 16° + 5°). We calculated the flux- 
weighted average by summing the fluxes of the individual stars 
contained in the aperture for each channel and calculating Q and 
U (and subsequently p and ff) from these total values. This cor- 
responds to a flux-weighted average over Q and U. 
The reason why this value is much lower than the polarization 
generally found in this area is the contribution of IRS IW with 
its high flux and an intrinsic polarization which is almost per- 
pendicular to that of the sources in the vicinity. 
This additionally supports our findings of higher polarization de- 
grees toward the eastern edge of the FOV: for the total polariza- 
tion to be on the order of 4% (including IRS 1 W), the surround- 
ing sources must have a significantly higher polarization degree. 
We compared our findings to several older NACO data-sets 
with the same FOV. Using aperture photometry on the IRS 1 
sources (except IRS 1 W) and the northern IRS 16 sources except 
IRS 16NW (to avoid problems with saturation), we found aver- 
age offsets in both polarization degree and angle in the order of 
10-15° and 2-3%, with higher polarization degrees and steeper 
polarization angles for the IRS 1 sources. This again confirmed 
the trends we find in our main data-set. 

Very few H-band polarization measurements are available for a 
comparison. The most recent surve y with an ap e rture not ex- 
ceeding our FOV was conducted by iBailev et all (1 19841) in the 
J-, H- and K-band, who present polarization parameters for two 
sources in our FOV, IRS 1 and IRS 16 (treating these complexes 
as a single source each, using a 3.0" aperture). That study mea- 
sured 9.9 + 0.6% at 20° ± 1° for IRS 1 respectively 10.2 ± 0.2% 
at 8° + 1° for IRS 16 in the H-band. Calculating a flux weighted 
average in a 3.0" aperture around the IRS 1 resp. IRS 16 sources 
based on our data yields values of 10.9 ± 0.5% at 15° ± 10° for 
IRS I and 7.2 ±0.5% at 21 ° + 10° for IRS 16. This agrees weU for 
IRS 1, while the polarization degree found for IRS 16 deviates 
considerably. But it has to be considered that both sources were 
only very poorly resolved at the time of that study, not all IRS 1 6 
sources are contained in our FOV, and we only used sources with 
reliably measured polarization for the comparison. 
Our measured polarization degrees and angles are also com- 
patible to the l arger- scale polarization maps presented by 
iNishivama et al.l (|2009) in the H- and the Ks-band. The authors 
find polarization angles of ~20° in the area around our own FOV, 
which itself is not covered in that study. 



3.4. Relation between H- and Ks-band polarization 

By comparing the positions of the sources detected in each 
datasets, we found 133 sources with reliable polarization pa- 
rameters common to the H- and the Ks-band data. The miss- 
ing sources are mostly found outside the other data-set's FOV, 
in addition to a small number of very fast moving sources (like 
the S stars), which are difficult to identify owing to the time of 
2.5 years between the H- and Ks-band observations. In addition, 
photometric errors tend to be larger in the westernmost region 
of the FOV because of the lack of suitable close bright stars for 
PSF determination. This leads to less sources with reliable polar- 
ization parameters there. The common sources are used for the 
following source-by-source comparison. 

We find that the polarization angles measured in the H- and Ks- 
band agree well in the eastern part of the FOV, while there ap- 
pears to be an offset in the center and the western region. FigfT2l 
shows the difference in polarization angle and the relation of H- 
and Ks-band polarization degrees. Oh-Oks can be fitted well with 



Table 1. Results of Gaussian fits to polarization parameter histograms 
for the complete dataset resp. sub-datasets separated based on p^s resp. 
position along the East- West-axis. Polarization degrees given in %, an- 
gles given in degrees. 
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Fig. 13. Sources detected in H- and Ks-band, separated by Ks-band polarization. White columns: pxs < 6%. Red columns: p/c, > 6%. Upper/lower 
left frame: Ks-band polarization degree/angle, Upper/lower middle frame: H-band polarization degree/angle, upper right frame: relation of H- and 
Ks-band polarization degree (logarithmic plot). Lower right frame: difference between H- and Ks-band polarization angle. Blue lines represent 
Gaussian fits to the histograms. 



a Gaussian distribution and shows a peak at 4° + 8° (FWHM of 
20°). Considering the width of the peak, this offset is not sig- 
nificantly different from zero. -2^ can be fitted quite well with 
a log-normal distribution, peaking at 1.9 ± 0.4 (FWHM of 0.9). 
For a complete list of the fitting results, see TablT] (which also 
contains the values referred to in the paragraphs below). 
Assuming the two peaks we find for both the polarization degree 
and angle in the Ks-band are real, we separated the stars with re- 
liable polarization parameters in both H- and Ks-band into two 
samples: stars with p^s < 6% (pK") and with pKs > 6% (pK^). 
FigfT3] shows histograms of the different polarization parameters 
of the pK" and pK^ sources in the two bands, namely the polar- 
ization degrees and angles, Oh-Qks and -^. All histograms were 
fitted with a Gaussian, and although these fits are poor in several 
cases, the fitted peaks show at least the trends present in the data. 
We find that the pK+ sources show systematically lower polar- 
ization angles than the pK" sources in the Ks-band (peak offset 
of 18°). The peaks fitted here match the ones fitted to the com- 
plete dataset (see 33.11 1. A similar, yet smaller offset exists in the 
H-band (7°, well within the uncertainties, the fitted peaks also 
correspond to those determined in 33.2b . Accordingly, Oh - Oks 
also shows an offset of 9°, but that value is relatively close to 
zero for both sub-datasets considering the large FWHM of the 
peaks. 

Looking at the polarization degrees, we find similar offsets: for 
the Ks-band, the pK^ peak is found at a polarization degree 
which is higher by a factor of 1 .6 than where the pK" peak is fit- 
ted. The relative difference found in the H-band is smaller, with 
the pK^ peak found at a polarization degree which is higher by 
a factor of 1 .3 than that of the pK" peak. This manifests itself in 
the -^ histogram, which we find peaked at 1.6 for pK^, and at 
2.0 for pK-. 

Because we mostly find the higher polarized pK^ sources in the 
eastern part of our FOV (in the general area of the Northern 



Arm), we also divided our sources into two samples based on 
their position along the East- West-axis: pK-E (sources more than 
4.1" east of Sgr A*) and pK-W (sources less than 4.1" east of Sgr 
A*). FiglT4]shows histograms of the polarization parameters for 
both sub-datasets. We find results that are very similar to those 
found for a separation based on pKs, with practically identical 
peaks and offsets (pK-E corresponding to pK^, pK-W to pK ). 
How can we explain these findings? For interstellar polarization 
in general, the H-band polarization degrees are expected to be 
significantly higher than the Ks-band values, while the angles in 
both bands should be the same within the uncertainties. This is 
expected from the Serkowski l aw and the power law relation pre- 
sented bv lMartin et al.l (Il990l) . According to the semi-empirical 
Serkowski law, the polarization at a given wavelength in relation 
to the polarization maximum depends on the wavelength where 
that maximum occurs 



p{X) 



— exp 



-K X In^ 



(^max \ 



(6) 



with K = 0.01 + 1.66A,„ajc (IWhittet et al.Lll992h . It appears, how- 
ever, that polarization in the NIR, specifica lly in the J-, H- an d 
K-band, only very weakly depends on A„,a^ (iMartin et al.Lll990h . 
Keeping this in mind and considering the availability of only two 
data-points for each source and the large FWHM of the Gaussian 
fits to ■^, we can only give rough estimates here. We find that 
the peak fitted to the complete dataset (1.9 ± 0.9) agrees best 
with /l,„flf ~ 0.7/im. Of the two fitted peaks for the sub-datasets, 
the value of 2.0 + 0.7 would agree with A,„ax ~ 0.1 jjm as well, 
while the peak at 1.6 + 0.7 points to a Amax which is eithe r much 
smaller (~ 0.25yum) or larger (~ 1 .25fim). iBailev et al.l (11984 
give a value of A„ax x 0.80/im (using J, H-and K-band data), 
which approximately matches the first two values we find here. 
The authors state that the agreement with the semi-empirical law 
is only rough. It also has to be considered that this law was estab- 
lished for sources with only very weak polarization in the NIR, 
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Fig. 14. Sources detected in the H- and Ks-band, separated by position along the East-West-axis. White columns: sources less than 4.1" east of 
Sgr A*. Red columns: sources more than 4.1" east of Sgr A*. Upper/lower left frame: Ks-band polarization degree/angle, upper/lower middle 
frame: H-band polarization degree/angle, Upper right frame: relation of H- and Ks-band polarization degree (logarithmic plot). Lower right frame: 
difference between H- and Ks-band polarization angle. Blue lines represent Gaussian fits to the histograms. 



and consequently does not describe observations in that wave- 
length regime very well. All this leads us to conclude that our 
data clearly are not sufficient to give a reliable estimate for this 
parameter. 
Using the power-law relation proposed bv lMartin et al.l d 19901) 



Ph 

PKs 






(7) 



we obtain a power-law index of o- = 2.4 ± 1 .7 for -^ = 1 .9 + 0.9, 
while the two sub-dataset peaks lead to ff = 1.7+1.6 resp. 
a - 2.5 ± 1.3. These values agree with the range of 1.5- 
2.0 given bv lMartin et al.l ([1990), although the uncertainties are 
quite large. It has to be stressed that our values have been ob- 
tained on a relatively narrow region and that a study of a much 
wider region is necessary before any reliable conclusions can be 
drawn on this matter 

FigfTSJ shows binned plots of the H- and Ks-band polarization 
parameters. Only the common reliable sources were used for 
this comparison. We averaged the polarization parameters for all 
sources contained in 1.3" wide bins along the East- West axis. 
We observe the same large-scale trends found based on the his- 
tograms, and the first impression points toward the same effects 
being present in both Ks- and H-band for both polarization de- 
gree and angle. Looking at the plots of -^ and 9h-0ks, however, 
the differences between the two bands become apparent: as seen 
in Figflll -^ is found to be around 2-2.2 in the center and more 
toward 1.5-1.6 east of Sgr A*. A similar trend is visible toward 
the western edge, but there are few sources there and there may 
be an influence of edge effects. The Oh - Oks plot also confirms 
the trends found earlier: we see an offset of ~10° in the center 
and the west and one of ~0° toward the east. 
These results raise the question for the cause of the observed de- 
viations of the polarization parameters over the FOV, assuming 
those deviations are indeed real and not some sort of instrumen- 



tal effect we did not account for. The comparison of our main 
Ks-band dataset (2) with dataset 4 shows that while a similar 
increase in polarization degree exists there as well, no accompa- 
nying shift in polarization angle is found. If this was indeed an 
instrumental effect, one would expect the same pattern in both 
parameters. Fig |A.3l shows a residue of the dither pattern in both 
Q and U, but these variations are far too small to lead to a sig- 
nificant deviation over the FOV as we observe it here. 
It has to be considered that the significance of the trends we 
find stays below 3cr, based on the systematic uncertainties de- 
termined in 33.11 

This question cannot be solved completely without additional 
(calibration) observations, but for the purposes of this study, we 
assume that what we observe is indeed a real effect. As possible 
explanations, two basic mechanisms come to mind here: 



Variable LOS extinction The extinction toward the central par- 
sec is known to be "patchy" (see Schodel et al., 2010b), which 
in turn indicates different dust column densities and/or dust pa- 
rameters along individual lines-of-sight toward different regions 
in the FOV. This could lead to differences in the polarization 
measured at different locations. But the situation is even more 
complex. Not only the densities are important, but also possible 
different alignment in individual dust clouds (so passing through 
an additional cloud on one LOS compared to another could even 
lower the total polarization degree for that LOS). Another prob- 
lem is that we find a smaller effect in the H-band (and thus the 
impact on Oh-Oks and ■^). This would require significantly dif- 
ferent average dust parameters from one LOS to the other (which 
in turn requires even more dramatic changes for a considerable 
percentage of the individual dust clouds along the LOS). Such a 
configuration is possible, but the highly specific arrangement re- 
quired to produce a pattern as we observe it seems very unlikely. 
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Fig. 15. H- and Ks-band polarization parameters, averaged over 1.3" bins along the East-West axis. Upper left: polarization degrees, red: Ks-band, 
green: H-band. Upper right: polarization angles, red: Ks-band, green: H-band. Lower left: ■^^. Lower right: 9h - 9ks- 



Local influences The area where this effect occurs coincides 
with the position of a Icnown local feature, the Northern Arm of 
the Minispiral. This feature is clearly visible in the L-band (see 
Figlll, faint in Ks and not detectable in the H-band. The light 
from the stellar sources itself, which passes through the stream 
of aligned grains in the Northern Arm, and scattered and/or emit- 
ted light from these grains themselves may contribute to some 
extent to the polarization measured in this region. This would 
have a stronger impact in the Ks-band compared to the H-band, 
because of the sizes and temperatures of the involved grains, this 
far matching our findings. The question however remains how 
substantial such a contribution could be. 

Under conditions as they are found in the filaments, grain align- 
ment by the Davis-Greenstein mechanism would be almost per- 
fect, especially because of the str ong magnetic fields (lower limit 
of ~2 mG close to IRS IW, see lAitken etalL 1998). By coin- 
cidence, the local grain alignment angle in the Northe rn Arm 
match es that measured for the LOS polarization ( Aitken et al.L 
Il998l) . This means that local dichroic extinction would increase 
the observed total polarization degree, while scattering/emission 
would decrease it (as is indeed found for IRS 1 W, see below). 
While the trends in polarization degree in the 2009 data may 
be explained this way, the measured polarization angles do not 
show the expected behavior: if local and LOS dichroic extinc- 
tion produce alignment in the same direction, the polarization 
angles should follow the magnetic field lines (see Fig|5]l. This 



is the case between 5" and 1" west of Sgr A*, but farther west, 
the angles are ~ 20° steeper than what would be expected. In 
the H-band, the angles are less steep in that region, providing a 
better match to the field lines. The preliminary Ks-band results 
(March 2011) also seem to agree well with the field. Farther to 
the south, however, the higher polarization degrees in the 2007 
Ks-band data-set occur despite a complicated field structure that 
does not offer an easy explanation for this effect. 
It is questionable whether local dichroic extinction can account 
for these effects, or if other contributions from scattering or ex- 
tended emission play a role here. The region observed here is 
quite complex, with a complicated field structure and gas/dust 
streams, and in addition, the exact position of the observed 
sources along the LOS axis is not clear. It is therefore difficult 
to determine what effect even a given dust distribution with a 
known temperature and composition would have on a certain 
source, but not even these parameters are known exactly. 
We therefore conclude that local effects, such as local dichroic 
extinction, extended emission by and scattering on dust grains 
are the most likely cause for this effect, although we cannot de- 
termine the likely extent of the individual contribution from the 
available data on dust parameters and configurations along the 
LOS and in the GC itself. 

This immediately leads to the next question: are these trends mir- 
rored by the behavior of the extinction measured in the FOV? 
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Fig. 16. Polarization efficiency in Ks-band (upper frames) and H-band (lower frames), compared to Ks-band extinction, plotted as point density, 
with typical error represented in upper left corner pK^ and pK" sources shown separately in left resp. right frames for both bands. Green lines 
represent the fitted power-law relation. 



3.5. Correlation with extinction 

We compared the H- and Ks-band polari zation values to the ex- 
tinction map presented bv lSchodelet al.l (|2010b). FigfTSl shows 
the polarization efficiency j^ for both bands plotted against the 
Ks-band extinction A^s taken from the extinction map at the lo- 
cation of each source. As it turns out, almost the same distri- 
bution of Aks is found for the pK^ and the pK" sources (see 
FigfTTli. which in turn leads to an offset between the two sub- 
datasets in polarization efficiency. We therefore plot them sep- 
arately in Fig[T6] In both bands, the distributions can be fitted 
with a power law: 






ccA' 



Ks' 



(8) 



with Pks.- - -0.4 ± 0.4 resp. fin.- = -0.6 + 0.6 in the H-band 
for the pK" sources. Despite the large errors which stem from 
consider able scatter of the param eters, this matches the relation 
found by iGerakines et al. (Il995h for th e Taurus dark cloud and 
also the results of IWhittet et ail (l2008h for Ophiuchus (both in 
the Ks-band). This is not directly comparable, because the GC 
is obscured by more than one dust cloud wi th possibly d iffer- 
ent dust alignments. A more general study bv lJonesI (Il989l) . ex- 
amining a great number of sources covering a range in optical 



depth of about a factor of 100, finds a power law relation with 

/3 0.25. In this study, the author proposed a model where the 

magnetic field along the LOS consists of a constant and a ran- 
dom component (see also Heilis, 1987), thus leading to different 
grain alignment in each section along the LOS. This reproduces 




Fig. 17. Ks-band extinction for sources with pxs < 6% (open columns) 
resp. Pks > 6% (filled columns). 
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Fig. 18. Map of the intrinsic Ks-band polarization of the extended 
source IRS IW. The arrows indicate the proper motions of the central 
source, the motion of the Northern Arm material and the motion of both 
relative to each other. 



the findings in that study quite well, and it is also consistent with 
our own results within the uncertainties. 

For the smaller number of pK^ sources, we find similar power 
law indices: ySA-.s,+ - -0.5+0.7 resp. /?//.+ = -0.4+ 1.5. Compared 
to the pH/K" values, we find a significant offset in polarization 
efficiency, while the underlying power law appears to be very 
similar. This might indicate that the additional polarization is in- 
deed caused by a local contribution, likely of Northern Arm ma- 
terial. To produce this deviation along the LOS, a very specific 
and therefore unlikely dust configuration would be required. 

3.6. Examining the extended sources 

Our polarimetric data cover two bright extended sources in the 
central parsec, IRS IW and IRS 21. The former is contained in 
both our H- and Ks-band data (datasets 1 and 2), while the latter 
is covered by several rotated Wollaston datasets of poorer qual- 
ity (datasets 3-15). IRS 1 W shows a clear horseshoe shape as ex- 
pected for a bow-shock source in high-quality Ks-band images, 
while IRS 21 does not. Owing to its high apparent polarization, 
IRS 21 would be interesting as a polarimetric calibration source, 
but only if the polarization is not variable. To constrain the na- 
ture of this source even more, we also conducted a flux variabil- 
ity analysis on IRS 21 and other extended sources in the FOV of 
our polarimetric data and several NACO imaging datasets, taken 
between 2002 and 2009 in the H-, Ks- and L-band. 

3.6.1. IRS1W 

Source morphology IRS IW shows the characteristic 
horse shoe-shape of a bow-shock source (see e.g. ' Tanner et al.L 
l2005h . This shape can already be made out in the raw images, 
but it becomes even more apparent in a Lucy-Richardson 
deconvolved image (using a PSF obtained from bright IRS 16 
sources, see FigfTSl). The observed shape agrees very well with 
the re lative ve locities of the streaming material of the Northern 
Arm (iLacy et al., 1991 ) at the location of the sour ce and the 
proper motion of IRS 1 W itself dSchodel et al.Ll2009l , velocities 
plotted in FigfTSTl. 
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Fig. 19. Map of the intrinsic H-band polarization of the extended 
source IRS IW (see FigfTSll. 



Spatially resolved polarimetry We measured the total polariza- 
tion of IRS IW as (1.8 + 0.5) % at (-37 + 5)° East-of-North 
after application of the M iiller matrix to account for instrumen- 
tal polarization. lOtt et al.l (|T999) provide values of (4.6 + 2.5) 
% at (-85 + 8)° East-of-North. It has to be considered that in- 
strumental effects take place on the same order of magnitude as 
the measured polarization degree, which may explain the large 
offset in polarization angle compared to the older values where 
these effects were not compensated. We attribute our lower to- 
tal polarization degree to the better angular resolution and thus 
less influence from neighboring sources. The values provided by 
lEckart et alJ d 19951) are clearly influenced by neighboring stellar 
sources, with (3.0 + 1.0) % at (10 ± 7)°. 

If we assume that the foreground polarization for this source is 
the same as for the surrounding objects and apply a depolariza- 
tion matrix with the parameters for p = 7.6%, 6 - 9.2°, we find 
a total intrinsic polarization of (7.8 ± 0.5) % at (-75 ± 5)°. 
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Fig. 20. Ks-band (red circles), H-band (green squares) and L-band 
(blue open crosses) lightcurves of IRS IW. Shown here is the flux of 
the source normalized to the flux of non-variable reference sources in 
the respective band. H- and L-band values were multiplied by a constant 
factor to match the average Ks-band values. 
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Fig. 21. IRS 21 observed on 2004-08-30, right: image before Lucy-Richardson deconvolution, left: image after deconvolution (logarithmic gray- 
scale, contours trace 10, 20, 30, 50, 60, 70, 80, 90, 95% of maximum flux in the left frame, resp. 40, 50, 60, 70, 80, 90, 95, 99% in the right 
frame. 
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Fig. 22. IRS 21 observed on 2005-05-14, right: image before Lucy-Richardson deconvolution, left: image after deconvolution (logarithmic gray- 
scale, contours trace 10, 20, 30, 50, 60, 70, 80, 90, 95% of maximum flux in the left frame, resp. 40, 50, 60, 70, 80, 90, 95, 99% in the right 
frame. 



In the H-band we measured the total polarization of IRS IW as 
(5.2 + 0.5) % at (12 + 5)° East-of-North. The polarization an- 
gle appears typical for a stellar source affected by foreground 
polarization, but the polarization degree is much lower than the 
~12% found for stellar sources in the vicinity. We applied a de- 
polarization matrix with p — YL%, 0=15° and found a resulting 
intrinsic polarization of (6.9 ± 0.5) % at (-73 + 5)°. The angle 
agrees very well with the Ks-band polarization angle, suggesting 
that the same process is responsible. The lower intrinsic polar- 
ization degree points to the lower influence of the extended dust 
component compared to that of the central source at this wave- 
length. 

We measured the polarization of individual regions of the source 
in the deconvolved images and subtracted the foreground po- 
larization using a Miiller matrix as described in ^2.51 The results 
are shown in FiglTS] We find polarization degrees of about 10-20 



% with very similar polarization angles for regions with signifi- 
cant flux. Apparently, the polarization degree is lower by a factor 
of up to ~ 3 around the apex compared to the tails. 
We applied the same technique to our H-band dataset, and FigfT9l 
shows the result for the immediate area around IRS 1 W. We find 
a polarization pattern comparable to that in the Ks-band, with 
polarization degrees of about 10-20 % with very similar polar- 
ization angles for regions with significant flux and much lower 
polarization degrees in the central region. There is no signifi- 
cantly lower polarization degree toward the apex, as found in the 
Ks-band, but it has to be noted that the horseshoe shape is much 
less pronounced in the H-band anyway. The angles agree with 
those determined from the Ks-band data. 

The total intrinsic polarization and the spatially resolved pattern 
of IRS 1 W can be explained from the combination of the motion 
of the source itself, the streaming velocity of the Northern Arm 
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Fig. 23. Map of the intrinsic Ks-band polarization of the extended 
source IRS 21 (see FigfTSll. 
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Fig. 24. Ks-band (red circles), H-band (green squares) and L-band 
(blue open crosses) lightcurves of IRS 21. Shown here is the flux of 
the source normalized to the flux of non-variable reference sources in 
the respective band. H- and L-band values were multiplied by a constant 
factor to match the average Ks-band values. 



and the magn etic fields present in that structure. 
lAitken et al.l (il998) mapped the polarization of the Northern 
Arm at 12.5 /im and inferred the magnetic field orientation 
by assuming that the polarization was produced by emission 
from magnetically aligned elongated dust grains. The magnetic 
field at the location of IRS IW is perpendicular to the po- 
larizati on angles we find he re. The projected velocity of the 
source (ISchodel et al.L l2009l) is parallel to the field lines and 
also parallel to the streaming motion in the Northern Arm it- 
self (Lacvetal., 1991). This leads to the field lines following 
the morphology of the shock around the source. At the apex this 
causes a weakening of the field, while it is compressed in the 
tails. This in turn leads to a weaker resp. stronger grain align- 
ment in the apex resp. the tails. 

Gr ain temperatures in the Northern Arm reach up to ~200-300 
K dSmith et alill990l:lGezariLll992h . and this is sufficient to ex- 
plain the observed 12.5 fim emission. But for significant emis- 
sion in the H- and Ks-band, much higher grain temperatures of 
~1000 K would be required. This raises the question if the ex- 
tended emission in the bow-shock is emission from or scatter- 
ing on the aligned grains. One possibility is that the tempera- 
ture in the shock i s inde ed high enough for significant emis- 
sion. iGeballeietal] (l2004l) suggest that these higher temperatures 
could be reached by very small grains (~0. 001-0. 01 fim, which 
is by a factor of 10 smaller than typical grains expected here) if 
they are heated by occasional high-energy phot ons or stochas- 
tic col lisions with high-energy electrons or ions. lMoultaka et al] 
(l2004l) find that the spectrum measured for IRS IW (including 
the bow-shock) matches a 900 K blackbody. At this higher tem- 
perature, emission from heated aligned dust grains should con- 
tribute to the Ks-band and H-band polarization. 
This leaves the question of the survival of these very small 
grains in a bow-shock environment and at this temperature. In 
addition, these grains would have to be aligned to produce the 
observed polarization. But the very same processes that would 
heat the small grains to high temperatures would also random- 
ize any previous alignment unless the alignment mechanism is 
much faster than the randomization. This may be the case here, 
because strong frozen-in and compressed magnetic fields pro- 
vide an even stronger and faster alignment than in other regions 



of the Northern Arm. This might also add to the lower polariza- 
tion at the apex due to turbulence, which could lead to a par- 
tial randomization of grain alignment, before the stronger field 
and uniform streaming motion in the tails increase the alignment 
again. 

In addition to emission, scattered light from the central source 
enclosed in a dusty envelope could contribute to the observed po- 
larization. This process is known to occur in planetary nebulae 
and d usty young stellar obje cts (YSOs) (e.g. Lowe & GledhilJ, 
2007| iLucas & Rociii Il998h . In the ideal case of pure Mie- 
scattering on spherical dust grains, most of the light would be 
scattered forward, but a significant portion is scattered perpen- 
dicular to the incident direction. This latter part is linearly po- 
larized, with the polarization vector in the plane of the sky and 
perpendicular to the original propagation direction of the light. If 
a source like this is viewed face-on, no total intrinsic polarization 
is detected, because the polarizations of the regions surrounding 
the central source cancel each other out. If the source is inclined, 
any total polarization angle can be produced. Spatially resolved 
polarimetric measurements of these sources show a characteris- 
tic centrosymmetric pattern of the polarization vectors, however. 



Table 2. Polarization parameters measured for IRS 21. The average 
over all observations resp. the standard deviation is given in the last 
line. 





date 


p[%] 


dp[%] 
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2007-04-01 


10.4 
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Fig. 25. Ks-band (red circles), H-band (green squares) and L-band 
(blue open crosses) lightcurves of IRS lOW. Shown here is the flux of 
the source normalized to the flux of non-variable reference sources in 
the respective band. H- and L-band values were multiplied by a constant 
factor to match the average Ks-band values. 



Clearly, this is not the kind of pattern we find here. 
But what if the grains are not spherical, but elongated and 
aligned as is the case here at least for a significant part of the 
grain population? Lucas & Roche ( 1998) find patterns of aligned 
polarization vectors similar to the ones we find here in the cen- 
tral regions of a minority of the sources examined in their study. 
They claim that these patterns cannot be explained by scattering 
on spherical grains, but that aligned elongated grains must play 
a role there. Whitney & Wolff (2002) modeled scattering and 
dichroic extinction for non-spherical dust grains, finding a vari- 
ety of polarization patterns for different input values for optical 
depth, degree of grain alignment, and inclination of the source. 
In general, at low optical depths they find polarization vectors 
perpendicular to the axis of grain alignment (i.e. the orientation 
of the angular momentum vector of the spinning grains), while 
high optical depth leads to a predominance of dichroic extinction 
and thus to polarization vectors parallel to the axis of grain align- 
ment. This study is focussed on spherical dust configurations and 
disk-like structures, so the results are not directly applicable to a 
bow-shock source. 

In the light of these results, we consider it likely that scattering 
on elongated grains contributes to the observed patterns in the 
H- and Ks-band data. Both scattering and emission should pro- 
duce polarization at the same angle. Without thoroughly mod- 
eling the conditions in such a bow-shock environment, we can- 
not conclude which process is dominant, only that both probably 
contribute. 

Flux variability The total flux of IRS IW can be diflicult to 
determine, because the center of this source is saturated in 
many images. Unfortunately, the repairing algorithm included in 
StarFinder only works on point sources, because it assumes that 
the saturated source has a PSF similar to that of non-saturated 
sources in the vicinity. This is obviously not the case for an ex- 
tended source like IRS IW. We therefore only used data where 
the peak fluxes at the location of IRS 1 W stayed below the satu- 
ration threshold. Fig|20|shows the flux of IRS 1 W in relation to 
the reference flux in the H-, Ks- and L-band. The flux appears to 
be variable by ~30% in the Ks-band, while the H-band value is 
inconclusive. No clear periodicity is apparent, which might in- 



dicate either an eiTatic variability or a period on the order of or 
below the time resolution of our data. The flux variability ex- 
ceeds that found for IRS 16C (see Fig|4|, but this may in part 
be due to the difficult photometry on IRS 1 W. In the L-band, the 
flux rises to a maximum in 2005 and drops again by ~15 %. We 
consider this consistent with a constant flux within the uncertain- 
ties, especially since there is no correlated behavior over the H-, 
Ks- and L-band. In total, we cannot state whether or not IRS 1 W 
shows a source-intrinsic flux variability based on the available 
observations. 

3.6.2. IRS 21 

Source morphology Previo us studies show IRS 21 as a roughly 
circular, yet extended source jTanner et al.Ll2002ll200 5'). We find 
the same for our data, but after applying a Lucy-Richardson 
deconvolution, this changes significantly: Fig|2T| and |22] show 
IRS 21 before and after deconvolution for two datasets (2004- 
08-30 resp. 2005-05-14). There is no clear bow-shock shape vis- 
ible prior to deconvolution. After the LR process, there is still 
no clear horseshoe shape, but it is clear that the source is not 
circular in projection and thus most likely not spherical. The de- 
convolved images appear to contain a central source with a bow- 
shaped northern extension. It has to be noted that this shape is 
not this apparent in all our datasets, but the resolution that can 
be achieved on an extended feature like this critically depends 
on the data quality, especially in a very dusty environment like 
that of IRS 21. The LR algorithm also tends to 'suck up' flux of 
extended features into a central source (see AppendixlBl). For all 
but a few periods, however, the deconvolved images show an ex- 
tended structure at this location, with an east-western bar/bow to 
the north and a point-source-like feature to the south. We deem 
this shape to be consistent with a marginally resolved bow-shock 
type so urce. Our findings match those presented by Tanner et al] 
( |2004|) . who find a similar spatial structure and position angle of 
the bow-shock in the L-band. Furthermo re, the relative motion 
of the Northern Arm m aterial Lacv et all ( 1991) and the proper 
motions of IRS 21 (e.g. iTanner et all |2005[) would agree with a 
bow-shock in this direction (see projected velocities plotted in 
Figl^ 
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Fig. 26. Ks-band (red circles), H-band (green squares) and L-band 
(blue open crosses) lightcurves of IRS lOE. Shown here is the flux of 
the source normalized to the flux of non-variable reference sources in 
the respective band. H- and L-band values were multiplied by a constant 
factor to match the average Ks-band values. 
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Fig. 27. Ks-band (red circles), H-band (green squares) and L-band 
(blue open crosses) lightcurves of IRS 5. Shown here is the flux of the 
source normalized to the flux of non-variable reference sources in the 
respective band. H-band values were multiplied by a constant factor to 
match the average Ks-band values. 



From the images where the bow-shock shape was visible, we de- 
termined a projected distance between the peaks of the southern 
and the northern feature of ~ 400 + 100 AU (see FigEB- The 
projected distance was almost the same for all images, and the 
error given here is just a rough estimate. Considering a possi- 
ble inclination of the source, this constitutes a lower limit for 
the standoff dista nce between the ce ntral source and the apex 
of the bow-shock. Tan ner et aP (|2002') provide modeled standoff 
distances for candidate central objects, using the Ks-band radius 
they measure for IRS 2 1 of ~650 AU as an estimate for the stand- 
off distance. Using the radius makes this value an upper limit for 
the projected standoff distance. They find that this value agrees 
best with the following stellar types: an AGB type star, W-R, 
OfpeAVN9 or W-R WC9. Especially a Wolf-Rayet (W-R) type 
star would be able to produce a wide range of possible standoff 
distances (310-10100 AU). Our value agrees with these possibil- 
ities as well. 



Spatially resolved polarimetry IRS 21 is not covered by our 
main high-quality polarimetric dataset, so we measured its to- 
tal polarization in several other datasets with lower Strehl ratios 
(see Tabl3l datasets 3-15, the data quality is still sufficient for 
polarimetry on bright sources). The values we find for the indi- 
vidual datasets are shown in Tab|2] Considering the uncertain- 
ties, we conclude from these measurement that both the polar- 
ization degree and angle of IRS 21 can be treated as constant 
within our margin of error We find an average value of (9.1 
+ 0.2) % at (16.4 ± 0.3) deg rees, whi ch agrees very well with 
the polarization measured by Ott et al. ( 1999) and also with the 
polarization angle determined bv .Eckart et al.. (1995 ). The latter 
study mea sured a po larization degree which is ~60% higher than 
our value. lOtt et alJ (11999) explained this difference by different 
apertures used and that might indeed be the case. We did not find 
a clear indication of a variable polarization degree of IRS 21 in 
the available data, although several, low quality datasets led to 
some outlier values. Since this coincided with large variations in 
the polarization of other sources, we do not consider this a sig- 
nificant effect. 

After applying a depolarization matrix with polarization param- 
eters determined from the surrounding point sources (5% at 30°, 



determined on two sources close to IRS 21 and thus not as reli- 
able as the values for IRS 1 W), IRS 21 appears to have a total in- 
trinsic polarization of (6. 1 ± 0.5) % at (5 + 5) degrees. The polar- 
ization degree is slightly lower than that of IRS 1 W, while both 
sources have very different intrinsic polarization angles. Again, 
the angle we find is perpen dicular to the magne tic field orienta- 
tion in this region given bv I Aitken et alJ (Il998l) . Unfortunately 
no polarimetric H-band data covering IRS 21 is available. 
Fig|23]shows the polarization of individual regions. We find po- 
larization degrees of about 3-8 %, with less uniform polarization 
angles than those found for IRS IW in regions with significant 
flux. We also detect an increase instead of a decrease of the po- 
larization degree toward the apex. The polarization angle shows 
more variation compared to IRS 1 W. 

Again, the relation of the source motion and the local magnetic 
fields and streaming motions offer an explanation. The source 
moves almost perpendicular to the field, so the field lines are 
likely compressed in front of the shock and diluted toward the 
flanks. This leads to a higher polarization at the apex. If this field 
orientation is indeed preserved over the whole structure, it could 
explain the observed polarization pattern, because the polariza- 
tion angles would be perpendicular to the field lines. This would 
be the expected behavior for emission/scattering producing the 
polarization. 

For this grain alignment dichroic extinction would be expected 
to produce polarization angles perpendicular to what is observed 
here. This would reduce the overall intrinsic polarization, but it is 
questionable if the optical depth of the dust surrounding IRS 21 
is sufficient for a significant contribution of this process. Again, 
a detailed model of bow-shock polarization using measurements 
at different NIR wavelengths as input parameters might clarify 
the extent of the influence of these processes. For the purposes 
of this study, we assume that scattering and emission play the 
dominant role. 



Flux variability Fig|24l shows the flux of IRS 21 in relation to 
the reference flux. For the H- and the Ks-band data the flux 
significantly and steadily increases by about 50% over the 
observed period (June 2002 to May 2008). This corresponds 
to an increase in brightness of about 0.4 mag. That we find 
the same increase for both bands suggests that it is intrinsic 
to the source itself, because the H-band is dominated by the 
stellar component, while the Ks-band source is dominated by 
the dust enshrouding the stellar source. We find a lower (~20%) 
correlated flux increase in the L-band data as well. What could 
cause this increase in flux? 

From the available data we cannot determine whether this is 
a periodic variability or not because we only observe a steady 
increase in fl ux. This wo uld be possible for all candidates 
presented by iTanner et alj (20_02): both AGB-stars and Wolf 
Rayet stars may show such an increase in luminosity in this 
time-frame, especially if they are in a mass-losing phase. 



3.6.3. IRS 10 

IRS low shows a flux variability in all three bands, but there is 
no clear periodic behavior (see Fig|25] only values not affected 
by saturation are shown). The Ks- and H-band fluxes vary by 
about ~30%, while the L-band flux only shows about 10% vari- 
ability. We do not consider this a reliable detection of intrinsic 
variability, because no systematic trends are visible at the ob- 
servable timescale. By comparison, IRS lOE shows a strong and 
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Table 3. Details of thie observations used for this worlc. N is the number of exposures that were taken with a given detector integration time 
(DIT). NDIT denotes the number of integrations that were averaged onhne by the read-out electronics during the observation. The Strehl ratio was 
measured using the Strehl algorithm of the ESO eclipse package (see N. Devillard, "The eclipse software", The messenger No 87 - March 1997, 
publicly available at http://www.eso.org/pwjects/aot/eclipse/distrib/index.html given here is the average value over all images in each dataset. 



program 



date 



band mode 



channels N/channel NDIT DIT[sec] camera Strehl 



1 


073.B-0084(A) 


2004-06-12 


H 


2 


083.B-0031(A) 


2009-05-18 


Ks 


3 


179.B-0261(A) 


2007-04-01 


Ks 


4 


179.B-0261(A) 


2007-04-03 


Ks 


5 


179.B-0261(A) 


2007-04-04 


Ks 


6 


179.B-0261(A) 


2007-04-05 


Ks 


7 


179.B-0261(A) 


2007-04-06 


Ks 


8 


179.B-0261(D) 


2007-07-18 


Ks 


9 


179.B-0261(D) 


2007-07-19 


Ks 


10 


179.B-0261(D) 


2007-07-20 


Ks 


11 


179.B-0261(D) 


2007-07-20 


Ks 


12 


179.B-0261(D) 


2007-07-21 


Ks 


13 


179.B-0261(D) 


2007-07-23 


Ks 


14 


179.B-0261(D) 


2007-07-23 


Ks 


15 


179.B-0261(D) 


2007-07-24 


Ks 



Wollaston 0, 45 30 

Wollaston 0, 45 143 

Wollaston (rotated) 0, 45 15 

Wollaston (rotated) 0, 45 70 

Wollaston (rotated) 0, 45 23 

Wollaston (rotated) 0, 45 70 

Wollaston (rotated) 0,45 51 

Wollaston (rotated) 0, 45 130 

Wollaston (rotated) 0, 45 70 

Wollaston (rotated) 0, 45 70 

Wollaston (rotated) 0, 45 70 

Wollaston (rotated) 0,45 51 

Wollaston (rotated) 0, 45 124 

Wollaston (rotated) 0, 45 41 

Wollaston (rotated) 0, 45 30 



1 


30 


4 


10 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 


2 


15 



S13 


0.171 


S13 


0.272 


S13 


0.126 


S13 


0.217 


S13 


0.173 


S13 


0.103 


S13 


0.138 


S13 


0.188 


S13 


0.075 


S13 


0.084 


S13 


0.137 


S13 


0.137 


S13 


0.079 


S13 


0.023 


S13 


0.081 



clear periodic variability with the flux increasing by a factor of 
3 and a period of about two years (see Fig|26]l. The variability is 
correlated in the H-, Ks and L-band. This source has been clas - 
sified as a late-t ype Mira variable bv lTamura & Werner I (Il996h . 
while lOtt et al.l (Il999.) label it as a long-period variable. Our 
measurements agree with these classifications. 



3.6.4. IRS 5 

IRS 5 also shows a flux increase in the Ks-band between 2002 
and 2008 (see Fig|27] one value affected by saturation was taken 
out), although the increase is not as strong and clear as that found 
for IRS 21. The H-band data is much less clear for this source, 
and while it appears to be variable at that wavelength as well, the 
variability seems quite erratic. We therefore cannot conclude that 
this is an intrinsically variable source on the observed timescale. 



4. Conclusions 

We draw the following conclusions: 



1 . Like several previous polarimetric surveys of the central par- 
sec at much lower resolution, we find that the polarization 
follows the Galactic plane in general. This confirms that the 
polarization can indeed be attributed to aligned dust grains 
between the Galactic Center and the observer (foreground 
polarization). We find systematic and individual deviations 
from this general pattern, however Over our FOV we find 
that the polarization degrees and angles change toward the 
eastern region around IRS 1. This behavior occurs in the 
H- and the Ks-band data, with stronger relative deviations 
in the Ks-band. That can be expected if this is indeed a lo- 
cal effect produced by dust grains in the Northern Arm, be- 
cause typical temperatures and grain sizes in this region lead 
to a much more significant impact of radiation from this 
grain population at longer wavelengths. This is supported 
by a systematic change of -^ over the FOV, with values 
of 2.0+0.3 in the center and 1.6+0.3 in the Northern Arm 
region. The higher of these values leads to a polarization 
maximum at ~ O. T^um, which match es previous findings of 
■Ima^f ~ 0.8//m bv lBailev et al.l (I1984 . based on JHK obser- 



vations of a more extended region in the GC. This value rep- 
resents dust parameters, specifically typical grain sizes, and 
it can be expected that local contributions by distinct dust 
masses lead to a different value than the average along the 
LOS. Unfortunately, we cannot derive any more constraints 
for the local dust based only on these two bands, other than 
that the dust grains responsible for the local effects must be 
elongated and aligned (see Aitken et al., 1998), and that their 
temperature and size must be sufficient for their influence in 
the H/Ks-band. 

We find a correlation between the spatially variable extinc- 
tion toward the central parsec and the polarization efliciency, 
j^. This relation appears to be different for the two groups 
of sources we find in our FOV (sources only affected by 
LOS polarization respectively sources where local contribu- 
tions are significant): for the foreground-polarized sources, 
we find a power law relation with a similar power law index 
for the H- and the Ks-band (j6h = -0.6+0.6,j6^, = -0.4+0.4 
). This is comp atible with the measurements resp. the model 
by I Jones! (|1989), who assume a combination of a constant 
and a random component of the galactic magnetic field along 
the LOS to align the dust grains. For the population of higher 
polarized sources, we find a higher polarization efficiency 
in general, while the power law index for both bands re- 
sembles that of the sources with lower polarization (with 
/3h - -0.4 ± 0.7 andyS^j = -0.5 + 1.5, although it has to be 
cautioned that these fits are poor because of the small number 
of sources). The higher polarization efficiency points to the 
influence of local dust grains, which only contribute a small 
amount to the total extinction. But because these grains are 
aligned, they add a significant component to the total po- 
larization, because this alignment is not (partially) compen- 
sated by averaging over other dust populations with different 
alignment as is the case for every individual dust configura- 
tion along the LOS. These findings further support our pre- 
vious conclusion that the influence of local dust is quite sig- 
nificant not only in the MIR, but also in the H- and Ks-band. 
Intrinsic polarization takes place in several sources in the 
central parsec, not only at longer wavelengths as shown first 
by Knacke & Capps (|1977h . but also in the H- and Ks-band 
(e.g. Ott et al.K.199^ . With the Miiller calculus we can iso- 
late the intrinsic component for point sources and maps of 
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extended features. The resulting total intrinsic polarization 
angles for IRS 21 and IRS 1 W agree very well with what can 
be expected from source morphology and relative motion of 
gas/dust in the Northern Arm and the sources themselves. 
We find very similar intrinsic polarization degrees for both 
sources, with (7.8 + 0.5) % at (-75 + 5)° for IRS IW and 
(6.1 ± 0.5) % at (5 + 5) for IRS 21 (both Ks-band). Contrai-y 
to the wavelength dependency seen in the foreground polar- 
ization, we find a slightly lower H-band polarization degree 
for IRS IW compared to the Ks-band: (6.9 + 0.5) % at (- 
73 + 5)°. That the central source contributes a larger amount 
of flux in the H-band compared to the extended component, 
while the intrinsic polarization mostly stems from the bow- 
shock, explains this discrepancy. 

The spatially resolved intrinsic polarization patterns of both 
bright extended sources in our FOV do not show a cen- 
trosymmetric shape as would be expected for Mie scattering 
on spherical grains. Instead, the pattern found for IRS IW 
suggests that scattering on aligned elongated grains must 
play a significant role, possibly combined with thermal dust 
emission (in the Ks-band), because the grains in the bow- 
shocks may reach te mperatures on the order of 1000 K (see 
iMoultaka et al.L 12004,) . Both effects would produce polariza- 
tion parallel to the long axis of the grains and thus perpen- 
dicular to their angular momentum. It therefore appears that 
the angular momentum of the grains in the bow-shock sur- 
rounding IRS IW is aligned parallel to the relative motion 
of the source through the Northern Arm material. Taking 
the magnetic fields in the Northern Arm into account, this 
is most likely caused by the field lines being warped around 
the flanks of the bow-shock. This in turn compresses the field 
lines and thus increases the local field strength, which is then 
sufficient for fast magnetic alignment of the grains. 
For IRS 21 we find a different polarization pattern. It appears 
that the polarization vectors are mostly parallel to the rela- 
tive motion of the source through the surrounding material. 
Combined with the local field orientation, this again leads to 
the same dominant mechanism of scattering/emission as for 
IRS 1 W, if we assume that the field lines are compressed in 
front of the shock that is moving perpendicular to them. 

4. After deconvolution, IRS21 clearly does not show a circu- 
lar shape in the Ks-band. Instead, its shape is consistent with 
that of a poorly resolved bow-shock with its apex north of 
the central so urce. This i s consistent with the L-band shape 
presented by Tanner et al.l (|2004) and the expectations from 
the proper motions of the source in relation to the stream- 
ing motion of the Northern Arm. We find a standoff dis- 
tance of ~ 400 AU, this agrees with the upper limit given by 
iTanner et al. (2002). It is difficult to constrain the nature of 
the central source by this parameter alone, and most stellar 
types with strong win ds, such as Wolf-Ra yet stars, would be 
likely candidates (see lTanner et ali[2002h . 

5. We find a ffux variability for several bow-shock sources in 
the central parsec: IRS 21 shows a ffux increase of ~50% 
over 6 years in the H- and Ks-band (less, but still signif- 
icant in the Lp-band). Especially the H-band measurement 
indicates that this variability must be intrinsic to the central 
source and not only to the enve lope. The ca ndidates sug- 
gested for the central source by iTanner et al.l (12002) could 
all show this variability, especially a source that is currently 
in a mass-losing phase. 

While we also find a possible ffux variability for IRS 5, 
IRS IW, and IRS lOW in all three bands, no clear period- 
icity or trend is apparent for these sources at our time resolu- 



tion. We therefore cannot determine if these sources show a 
true intrinsic variability or not. We confirm the known peri- 
odicity of IRS lOE (~2 yrs) in the Ks-band, finding the same 
variability in the H- and Lp-band as well. 

Several effects appear to contribute to the observed polarization 
in the central parsec of the GC, in addition to the foreground 
polarization. The findings in our relatively narrow field-of-view 
alone indicate possible additional areas of interest: do the large- 
scale trends we find continue in other sections of the Northern 
Arm (and the other elements of the minispiral)? Do the other 
bow-shock sources, such as IRS lOW, 5, and 8 show similar 
polarization patterns, and what does that tell us about the lo- 
cal processes in the bow-shocks? Are other intrinsically polar- 
ized sources present in the field, and what type of objects would 
show this behavior? More observations with wider FOVs and in 
different bands (especially L-band, to directly examine extended 
dust emission e.g. in the Northern Arm) would be desirable, and 
the recently achieved direct polarimetric calibration of NACO 
offers the possibility to map more extended areas of the GC in 
polarization at high resolution. 
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Appendix A: Error estimation 

In general, several effects can contribute to the total flux un- 
certainty: the basic uncertainty in the counts measured for each 
pixel, the insufficiently determined PSF (resp. the fact that a PSF 
determined from several bright sources does not fit every source 
perfectly) and possible effects of insufficient calibration (e.g. er- 
rors in the sky or the flat-field). 
When using StarFinder, the formal flux error given by the al- 
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Fig. A.l. Stokes parameters Q and U (left frame) and resulting po- 
larization degree (right frame), measured for each pixel in a wiregrid 
flat-field. 



Fig. A.2. Stokes parameter Q measured at each pixel of a flat-field 
taken with (upper frames) and without (middle frames) the Wollaston 
prism. The lower frames show Q after dividing the Wollaston flat-field 
by the flat-field taken without the Wollaston. Left frames: H-band, right 
frames: Ks-band. 



gorithm represents the first of these three contributions. The al- 
gorithm assumes an independence of all the pixels resp. their 
measured counts, h owever. This is no t the case anymore for a 
deconvolved image. ISchodel I (1201 Oal) circumvented this prob- 
lem by introducing a modification factor for these formal errors 
(determined by simulations). Owing to the different instrument 
setup used here, we could not simply adopt the same factor 
The second contribution stems from a variation of the PSF over 
the FOV. Generally, the shape of the PSF changes with the rela- 
tive position of the source to the guide star (anisoplanasy), with 
longer projected distances leading to larger distortions. The ex- 
tent of this behavior is highly dependent on the data quality 
resp. the seeing. The coiTelation factor provided by StarFinder 
describes this effect, because it measures how well a source 
matches the PSF, but it is not clear how to translate this into 
a flux uncertainty. We therefore estimated this effect in another 
way: we extracted three different PSFs from the western, central, 
and eastern region of the FOV and used them for PSF fitting with 
StarFinder As it turned out, this does not only produce different 
fluxes in all four channels, but can also lead to a shift in polar- 
ization angle and/or degree (depending on the source). A similar 
effect can be seen when comparing fluxes and polarization pa- 
rameters measured in different sub-frames of the deconvolved 
image. 

We compared these measurements to aperture photometry ap- 
plied to two bright sources (IRS 16C, IRS IC). The results are 
shown in Figs jA.3IA.4l In general, we find that the values deter- 
mined by aperture photometry are matched better by the PSF- 
fitting values the closer the PSF reference sources are to the 
source in question. Thedeviationscanreachvaluesof up to 1.5% 
resp. 10° even for these bright sources. 

The reason for this behavior is revealed by a closer look at the 
individual PSFs: these show a polarization themselves, and that 
pattern changes over the FOV (see Fig jB.2i the polarization vec- 
tors are almost perpendicular in several locations, even in ar- 
eas with significant flux). This leads us to conclude that the 
fluxes determined with a "nearby" PSF are more accurate than 
those determined with a "distant" PSF, and we therefore apply a 
weighted average (based on the distance of the individual source 
to the average position of the PSF stars) to the photometry on the 
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Fig. A.3. Fluxes and polarization parameters of IRS 16C plotted against time of observation. Upper left: flux measured in four channels. Upper 
right: Q (black) and U (red). Lower left: polarization degree. Lower right: polarization angle. The plotted lines represent the polarization parameters 
obtained on the complete mosaic by aperture photometry (red), PSF fitting photometry using the IRS 16 stars (light blue) resp. the IRS 1 stars 
(magenta), resp. the final value after deconvolution assisted PSF fitting photometry (green). 



deconvolved sub-images. The results provide a better match to 
the aperture photometry (though it has to be cautioned that the 
latter is problematic in crowded regions and where large-scale 
extended emission is present, which is why this method is not 
applied to the image as a whole). The standard deviation of the 
fluxes over the sub-frames provides an estimate for the contribu- 
tion of this effect to the total uncertainty. 

Using these error contributions to estimate the final, total un- 
certainty is problematic: the sub-frames are all taken from the 
same mosaic and are thus not independent measurements, the 
formal errors given by StarFinder are unreliable and generally 
far too small for brighter sources, and possible calibration prob- 
lems (sky, flat-field) are not covered at all. 
We therefore used a different approach to determine the total 
flux errors: we applied the same photometric method that we 
used on the full mosaic to mosaics created from eight (Ks-band) 
resp. six (H-band) sub-datasets, selected based on the dither po- 
sition along the East- West axis. This led to sub-sets of about 15 
(Ks-band) resp. 5 (H-band) images per dither position, and the 
resulting mosaics have the same depth resp. signal-to-noise ra- 
tio everywhere. Because the original images can be regarded as 
independent measurements, so can the resulting mosaics, with 
the additional advantage that fainter stars still have a sufficient 
signal-to-noise ratio (as compared to individual images). The 
downside of this approach is that all sources in the eastern- and 
westernmost ~2" are contained in only one of these sub-mosaics. 
We then adopted the standard deviation of the fluxes measured 



in these sub-mosaics as the photometric error of each source. 
In several cases (mostly in the H-band) this method yielded ex- 
tremely small errors for weak sources (see FigO. This can hap- 
pen by coincidence if a source is only contained in a small num- 
ber of sub-mosaics. We consider these errors (and the photom- 
etry of these sources in general) unreliable, and the respective 
source will have a tendency not to pass the reliability criterion 
for the polarimetry. 

The aperture photometry applied to individual frames shown e.g. 
in Fig lA. 31 reveals two other effects: the measured flux varies 
over the frames, and in addition to the instrumental polariza- 
tion (most prominently the "bump" in p), small variations in Q 
and U (and subsequently p and 0) can be observed. These vari- 
ations coiTespond to the dither pattern used in the observation, 
and they probe the different areas of the Wollaston prism resp. 
the detector. These deviations are much smaller than the photo- 
metric errors, however, so we assume that we can safely neglect 
resp. average over t hem . 

IWitzel et al.l (1201 lb determined the system-intrinsic systematic 
uncertainty of polarimetric observations with NACO as 1% in 
polarization degree and 5° in polarization angle. For a source 
with 5% polarization (a typical value for GC sources in the Ks- 
band), this already corresponds to a relative error of 20%. We 
simulated the required accuracy of the photometry so that the 
final polarimetric errors stay below three times (resp. six times 
for the H-band) the systematic uncertainty. We tolerate higher 
uncertainties in the H-band because the higher polarization there 
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Fig. A.4. Fluxes and polarization parameters of IRS IC plotted against time of observation. Upper left: flux measured in four channels. Upper 
right: Q (black) and U (red). Lower left: polarization degree. Lower right: polarization angle. The plotted lines represent the polarization parameters 
obtained on the complete mosaic by aperture photometry (red), PSF fitting photometry using the IRS 16 stars (light blue) resp. the IRS 1 stars 
(magenta), resp. the final value after deconvolution assisted PSF fitting photometry (green). 



leads to a lower relative error for the same absolute uncertainty 
of the polarization degree. This led to error thresholds of 3% in 
the Ks-band and 6% in the H-band. Most of the sources brighter 
than 16 resp. 18 mag (Ks/H) have errors below these thresholds 
(see Figl3]l. This leads to the final error of the polarization pa- 
rameters being dominated by the statistical photometric errors. 



Appendix B: LR deconvolution: useful for extended 
structures? 

It is known that Lucy-Richardson deconvolution has a tendency 
to " suck up" surrou nding flux into a central point source (see 
e.g. ISchodel L l2010ah . This impedes or severely hinders the de- 
tection of faint sources close to bright ones, introduces astro- 
metric and photometric uncertainties and makes a reliable back- 
ground estimation impossible. But unlike linear deconvolution, 
this method allows flux measurements on extended structures via 
simple aperture photometry on the deconvolved images. Can it 
be used to measure the polarization of the extended sources in 
the GC? 

To investigate this, we produced images of an artificial bow- 
shock-like structure next to a point source. Several other point 
sources with different fluxes were added to simulate a crowded 
field as is found in the GC. For each filter (0°, 45°, 90° and 135°), 



we modified all fluxes by a factor calculated for generic polar- 
ization parameters of P = 4-5% and 6 - -10°-100°: 



fpoLrel =l+PX 003(9 - 5), 



(B.l) 



This was applied to the extended feature as well, with a constant 
polarization degree of 10% and the polarization angles parallel 
to the bow-shock structure. We then added an additional fore- 
ground polarization of 4% at -15° to the extended feature. The 
resulting map was convolved with a PSF extracted from the real 
images. We also added a weakly position-dependent background 
in the shape of a very flat two-dimensional Gaussian similar in 
level to that found in the actual observations. 
We applied the Lucy-Richardson deconvolution algorithm to 
these artificial images, followed by aperture photometry on the 
point sources and the extended structure, mapping the latter 
with two pixel apertures along its length. The structure itself 
was successfully recovered, although slightly widened. Fig jB.ll 
shows a comparison of the recovered polarization parameters 
in comparison to the input values. It appears that the only 
significant large deviation occurs for the source in the center 
of the extended feature: the recovered polarization degree is 
2.6 percentage points lower than the input value, while the 
polarization angle is 5° less. This appears to be an eff'ect of 
the algorithm "sucking" flux of the extended feature into the 
central source, thereby changing its polarization because of the 
significantly different polarization of the shell. The polarization 
degrees of several of the very faint point sources were also 
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Fig. B.l. Recovered polarization parameters for point sources (left frames) and extended feature (right frames). The green lines indicate a deviation 
of ±1.5% resp. 5° from a perfect recovery, indicated by the black line. 



recovered poorly, while at the same time their recovered polar- 
ization angles agreed much better with the input values. 
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Fig. B.2. Polarization pattern of the PSF used for the first step of the 
deconvolution-assisted photometry. Red: PSF determined from IRS 16 
sources. Green: PSF determined from IRS 1 sources. 



The values recovered for the bow-shock-feature itself show 
only small deviations compared to the input values, with a 
maximum difference in polarization degrees of 1.5-2% (smaller 
in most cases) and 5-7° in polarization angle. The deviations 
are even smaller for the brighter stars, with less than 0.5% and 
less than 1° in most cases. The fainter stars show only slightly 
larger deviations. This approaches the systematic uncertainties 
inherent to the instrument. We can therefore consider the 
Lucy-Richardson algorithm to be sufficiently accurate for our 
purposes of detecting extended features. The small variances 
found for the artificial "bow-shock" also indicate the level of 
accuracy that can be expected for the polarization parameters 
determined in ^3.61 

Fig lB. 21 shows the polarization pattern of two PSFs determined 
from different bright sources in the FOV (IRS 16 and IRS 1 
point sources). There are significant spatial variations, so 
these variations must occur (to some extent) in all sources 
in the non-deconvolved images. Because these sources are 
expected to be essentially point-like before their light enters the 
atmosphere, polarization patterns like this must be introduced 
by instrumental or atmospheric effects. This is another reason 
why deconvolution is an important step to recover polarization 
patterns of extended sources successfully: the intrinsic patterns 
would otherwise be masked by the strong spatial variations, 
with local values for the polarization degree of up to 15% in the 
PSF pattern (although these values are much lower in the central 
region with the highest flux, on the order of 1-2%). These effects 
are removed by the deconvolution process. 
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